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• ABSTRACT 

We study the 3-dimensional distribution of non virialised matter at z ~ 2 using high 
resolution spectra of QSO pairs and simulated spectra drawn from cosmological hy- 
drodynamical simulations. We have collected the largest sample of QSO pairs ever 
observed with UVES at the ESO-VLT, with angular separations between ~ 1 and 14 
' arcmin. The observed correlation functions of the transmitted flux in the H I Lyman a 

Q-f forest along and transverse to the lines of sight are in good agreement implying that 

q the distortions in redshift space due to peculiar velocities are small. The clustering 

S_h ■ signal is significant up to velocity separations of ~ 200 km s _1 , or about 3 hr 1 co- 

"^2 \ moving Mpc. The regions at lower overdensity (p/p < 6.5) are still clustered but on 

& ■ smaller scales (Aw < 100 km s _1 ). The observed and simulated correlation functions 

are compatible at the 3er level. A better concordance is obtained when only the low 
. t-H , overdensity regions arc selected for the analysis or when the effective optical depth of 

■ the simulated spectra is increased artificially, suggesting a deficiency of strong lines 

] in the simulated spectra. We found that also a lower value of the power-law index 

of the temperature-density relation for the Lyman a forest gas improves the agree- 
ment between observed and simulated results. If confirmed, this would be consistent 
with other observations favouring a late Hen reionization epoch (at z ~ 3). We re- 
mark the detection of a significant clustering signal in the cross correlation coefficient 
at a transverse velocity separation Awj^ ~ 500 km s -1 whose origin needs further 
investigation. 

Key words: intergalactic medium, quasars: absorption lines, cosmology: observa- 
tions, large-scale structure of Universe 
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1 INTRODUCTION 

The study of the evolution of cosmic structures is one of the 
major topic of present-day cosmology. 

It is now more than a decade that thanks 
to semi-analytical and hydrodynamical simulation 
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results (e.g., 


ICen et alJ 11994 IZhane et all 


19951: 


Hernauist ct al. 


119961: iMiralda-Escude et alJ 


19961; 


Bi & Davidson 


19971 iDave et alJ Il997l IZhane et alJ 


1997; 


Theuns. Leonard & Efstathioul Il998al iMachacek et alJ 


2000l) the majority of the absorption features observed in 



Lyman a forests of high redshift QSO spectra is identified 
with the fluctuations of the intermediate and low density 
intergalactic medium (IGM), arising naturally in the hier- 
archical process of structure formation. The physics of this 
highly ionized gas is simple, governed mainly b y the Hubble 
expan sion and the gravitational instability iRauch et alJ 
2005). If photoionization equilibrium is assumed, a mean 
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relation between tem perature and overdensity is obtained 
<Hui fe Gnedmlll997l) : 



po wer spectrum of the spatial distribution of dark mat- 
ter JCroft et alJll998l.ll999l.l2002l:lNusser fc Haehnerdll999l. 



T — To (p/p) 



(1) 



where Tb and 7 depend on the ionization history of the Uni- 
verse (for an early reionization, To = 10 4 K and 7 = 1.6 
are generally adopted) . The combination of photoionization 
equilibrium and eq. (1) leads to a power law relation between 
measure d Hi Lyman a optic al depth and gas overdensity 
(see e.g. IWeinberg et al.| [l999. for a review), which, adopt- 
ing the standard values for the parameters, can be written 
as: 



T"HI 



1.3 x 10^(1 + z) 9/2 (p/p) 



(2) 



A critical test of the nature of the Lyman a absorbers 
as proposed by simulations, comes from the determination 
of their sizes and spatial distribution. Correlations in the 
Lyman a forest were detected with a 4-5 a confidence by 
various authors at typical scales Av < 350 km s -1 observi ng 
at high r esolution individual l ines of sight ihu et al|l996l at 
z ~ 3.7; ICristiani et alJll997l at z ~ 3; iKim et aljhoOll at 
z ~ 2). This velocity range corresponds to scales < 2.5 h^ 1 
Mpc (considering peculiar velocities negligible). In this ob- 
servational approach, however, the 3-dimensional informa- 
tion is convoluted with distortions in the redshift space, due 
to peculiar motions and thermal broadening. 

Multiple lines of sight offer an invaluable alternative to 
address the spatial distribution of the absorbers, enabling 
a more direct interpretation of the observed correlations. 
The interesting range of separations lies about the Jeans 
scale of the photoionized IGM (~ 1 arcmin or ~ 1.4 h -1 
comoving Mpc at z ~ 2). At this scale there should be 
a transition from a smooth gas distribution, which pro- 
duces nearly identical absorption features in neighbouring 
lines of sight, to a correlated density distribution, where 
the correlation strength decr eases with incre asing separa- 
tion of the lines of sight (e.g., IViel et aLll2002T) . Indeed, the 
spectra of multiple images of len sed quasars with separa- 
tions of the order of a few arcsec JSmette et alJll99l Il995t 
llmpev et al1ll996t iRauch et all l200lft show nearly identi- 
cal Lyman a forests, implying that the absorbing objects 
have large sizes (> 50 h" 1 kpc). For pairs with larger 
separation the correlation between the absorption features 
becomes weaker. Studies based on the statistics of coin- 
cident and anti-coincident absorption line s provided evi- 
dence for dimensions of a f e w hundred k pc ( Bcchtold e^^| 
19941: iDinshaw et all Il99l Il995l Il997t iFang et alJ Il996t 
" Crotts fc Fandll998l: iD'Odorico et alJll998t iPetitiean et all 
19981: lLopez. Haeen fc Refiners! l200oi : lYoung et al.l l200lf 
Becker et al.l 120041) . The obtained large sizes were con- 
clusive to exclude models of the Lyman a forest as 
a population of p ressure confined small clouds (e. g., 
ISargent et alJll 98Cf) or clouds in d ark matter mini haloes 
(e. g.. lMiralda-Escude fc ReesHl993|) . 

The new theoretical scenario and the production of 
exceptional quality QSO spectra by high resolution spec- 
trographs at 10-m class telescopes, started the use of the 
Lyman a forest as a cosmological probe. From the study 
of absorption spectra along single lines of sight to dis- 
tant QSOs and the implementation of various algorithms 
to convert the observed flux into density it was possible, 
for example, to determine the shape and amplitude of the 



120001: IViel. Haehnelt fc Springein2004ah or, more straight- 
forwardly, of th e transmitted flux llMcDonald et alJ 120001 : 
IKim et al.ll2004h . 

The final goal of this paper is to investigate the dis- 
tribution properties of matter in the IGM applying the 
modern interpretation of the Lyman a forest. In particular, 
we computed the correlation properties of the transmitted 
flux along and across the lines of sight for a new excep- 
tional sample of high resolution spectra of close QSO pairs 
and group s at z eul ~ 2. Previous works following a s imilar 
approach jRollinde et alJl2003l : ICoppolani et al.l |2006l) used 
sets of low resolution spectra of QSO pairs with the main 
aim of deriving constraints on Oqa from the applicati on of 
the Alcock-Paczyhski test jAlcock fc Paczvnskilll979l see 
also Section 6). However, about 13(A#/(1 arcmin)) 2 QSO 
pairs with separation < A8 are needed to determine Qq a 
with a precision better than 10 per cent (lMcDonaldll2003fl . 
So, we will not attempt to determine f2oA in this study, al- 
though our plan is to gather a larger sample of QSO pairs 
in a few years. 

The present investigation is completed by the compari- 
son with a set of mock spectra drawn from an hydrodynam- 
ical simulation, reproducing several different realisations of 
our sample of QSO associations. We analysed how variations 
of the mean flux (corresponding to variations of the intensity 
of the ionizing background) and of the temperature-density 
relation of the gas affect the correlation functions. 

In this preliminary study, we do not take advantage of 
the presence, in our sample, of two groups of QSOs (one 
triplet and one sestet) to detect structures extending on 
large scales. Since new observing time have been allocated to 
our project at UVES-VLT, a forthcoming paper will be de- 
voted to the detailed analysis of coincidences and anticoinci- 
dences both of Lyman a and metal absorptions among mul- 
tiple lines of sight, on improved signal-to-noise ratio (S/N) 
spectra. 

The paper structure is as follows: in Section 2 we de- 
scribe the observed data sample, the reduction procedure 
and the simulated spectra. Sections 3, 4 and 5 are devoted 
to the computation of the auto and cross correlation func- 
tions of both the observed and simulated spectra and to their 
comparison. We discuss some interesting aspects of our re- 
sults and we explore the parameter space of the simulations 
in Section 6. The conclusions are drawn in Section 7. 

Throughout this paper we adopt flo m — 0.3, JIoa = 0.7, 
and h = H /72 km s" 1 Mpc" 1 . 



2 DATA SAMPLE 

2.1 The observed spectra 

The exploitation of the potential offered by multiple QSO 
lines of sight has been limited by the dearth of suitable 
groups of QSOs close and bright enough to permit high res- 
olution spectroscopy. Two major breakthroughs have dra- 
matically improved this situation: 

- the 2dF QSO Redshift Survey JCroom et alJ 120041. 
www.2dfquasar.org), whose complete spectroscopic cata- 
logue contains more than ~ 23000 QSOs in a single homo- 
geneous data base, which is approximately 50 times more 
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Figure 1. Relative positions, U magnitudes and redshifts of the QSOs composing the sestet (left plot) and the triplet (right plot) in our 
sample 



Table 1. Characteristics of the observed QSO spectra 





Object 


2 


M B 


Lyman a 
range 


S/N 
per pixel 


Pair A 


PA1 


2.645 


19.11 


2.094-2.585 


8-12 




PA2 


2.610 


19.84 


2.094-2.550 


3.5-6.5 


Triplet 


Tl 


2.041 


18.20 


1.633-1.991 


3-10 




T2 


2.05 


18.30 


1.592-1.999 


3.5-8.5 




T3 


2.053 


18.10 


1.665-2.002 


2.5-6 


Sestet 


SI 


1.907 


19.66 


1.665-1.859 


2-5.5 




S2 


2.387 


19.53 


1.858-2.331 


3-7.5 




S3 


2.102 


19.31 


1.633-2.051 


4-10 




Si 


1.849 


19.59 


1.575-1.802 


2-8 




S5 


2.121 


18.85 


1.633-2.069 


3.5-8 




S6 


2.068 


20.19 


1.592-2.017 


3-8.5 


Pair U 


UM680 


2.144 


18.60 


1.653-2.092 


6.5-17 




UM681 


2.122 


19.10 


1.634-2.070 


7-17 


Pair Q Q2343+12 


2.549 


17.00 


1.994-2.490 


13-23 




Q2344+12 


2.773 


17.50 


2.183-2.711 


12-18 



than the previous largest QSO survey to a similar depth 
(B<21); 

- the UVES spectrograph (Dekkcr ct al. 2000) at the Kueyen 
unit of the ESO VLT (Cerro Paranal, Chile) which has a re- 
markable efficiency especially in the extreme UV (close to 
the atmospheric cut-off). In this way, relatively low-redshift 
(z < 2.5) QSOs become accessible to high resolution obser- 
vations of the Lyman a forest down to faint apparent mag- 
nitudes and their surface density becomes high enough to 
provide several lines of sight for a tomography of the IGM. 

We searched the 2dF QSO database for the best groups 
with apparent magnitude B ^ 20 and z > 1.8. A triplet 
found in the Calan-Tololo QSO survey jMaza et al.lll993l 
Il995h was added to the sample. A great observational ef- 
fort was carried out to collect UVES spectra of the selected 
QSOs, which all have magnitudes fainter than B ~ 19, with 



the exception of the triplet for which B ~ 18. ESO allocated 
51 hours of observation to our project up to now, which al- 
lowed us to obtain acceptable S /N spectra of one pair (from 
now on called Pair A, with angular separation of ~ 1.3 '), 
one sestet and the triplet of QSOs (see Fig. . 

Reduction of the QSO spectra wa s conducted with the 
pipeline of the instrument (version 2.1, iBallester et all2000F) 
provided by ESO in the context of the data reduction pack- 
age MIDAS. In most cases, we could not apply the standard 
procedure due to the faintness of the spectra in the blue 
region and had to pre-filter the cosmic rays to make the 
optimal extraction work properly. Single extracted spectra 
were summed and rebinned and wavelengths were corrected 
to the vacuum-heliocentric reference frame. The final spec- 
tra have resolution R ~ 40000 in the Lyman a and C IV 
forest, while the S/N per pixel varies on average between 
S/N ~ 3 and 10 in the Lyman a forest and between 4 and 
15 in the Civ forest (see Table^for details). 

Continuum determination, in particular in the Lyman a 
forest region, is a very delicate step in the process of spectra 
reduction. Tentative procedures realised up to now to objec- 
tively determine the continuum position through automatic 
algorithms do not give satisfactory results. We adopted a 
manual subjective method based on the selection of the re- 
gions free from clear absorption that are successively fitted 
with a spline polynomial of 3rd degree. 

The limitations introduced by the uncertainty in the 
true continuum level and shape should be less important in 
the computation of the cross correlation function than in the 
case of single line of sight analysis. This because the undu- 
lations are uncorrelated between adjacent lines of sight, in 
part icular, if the two QSOs do not have too similar redshifts 
(see lViel et aljl2002l for a discussion in the case of power 
spectra). 

We added to the sample two more pairs from our 
archive: UM680/UM681 (Pair U, separated by ~ 1 arcmin) 
and Q2344+1228/Q2343+1232 (Pair Q, separated by ~ 5.57 
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Figure 2. Lyman cv-forests of two of the closest pairs in our 
sample. Pair A (upper two panels) separated by 1.3' and the pair 
at 1.06' formed by the QSOs S5 and S6 of the sestet (lower two 
panels) 



Table 2. List of the considered QSOs from the ESO Large Pro- 
gram (see text) 



Object 


Redshift 


M B 


Lyman a range 


HE 1341-1020 


2.134 


17.1 


1.64431- 


-2.08216 


Q 0122-380 


2.189 


16.7 


1.69072- 


-2.13625 


PKS 1448-232 


2.220 


16.96 


1.71688- 


-2.16674 


HE 0001-2340 


2.280 


16.7 


1.76750- 


-2.22575 


Q 0109-3518 


2.406 


16.44 


1.87381- 


-2.34966 


HE 2217-2818 


2.406 


16.0 


1.87381- 


-2.34966 


Q 0329-385 


2.423 


16.92 


1.88816- 


-2.36638 


HE 1158-1843 


2.453 


16.93 


1.91347- 


-2.39589 



arcmin). They were observed with LIVES at the same res- 
olution as the new data, and reduced following the same 
procedur e. More details on their properties are given in the 
paper bv lD'Odorico et al.l (I2002I) . 

The list of the QSOs and the main characteristics of the 
spectra are reported in TableQ The Lyman a forests of two 
of the closest pairs in our sample are shown in Fig. [5] 

The total sample is formed by 21 pairs uniformly dis- 
tributed between angular separations of ~ 1 and 14 ar- 
cmin, corresponding to comoving spatial separations be- 
tween ~ 1.4 and 21.6 Mpc. The median redhift of the 
Lyman a forest is z ~ 1.8. This is the largest sample of 
high resolution spectra of QSO pairs ever collected, unique 
both for the number density - we have 6 QSOs in a region 
of ~ 0.04 square degrees - and the variety of line of sight 
separations investigated. 

As a reference sample (see Sections 4 and 5.1), we con- 
sidered 8 LIVES QSO spectra obtained in the framework 
of the ESO Large Program 'The Cosmic Evolution of the 
IGM' dBereeron et al.ll2004i LP sample). The QSOs were 
selected in order to match as much as possible the redshift 



range of our sample (see Table |5J. The spectra have resolu- 
tion R ~ 45000 and S/N ratio ~ 50 in the Lyman a forest 
region. They were normalised to the continuum and pre- 
pared for the analysis following the same procedure adopted 
for the spectra of the pair sample. 

2.2 The simulated spectra 

In order to both assess the nature of the Lyman a forest 
inferred from simulations and to constrain the cosmological 
scenario of the same simulations, we compared the results 
obtained for our sample of observed QSO spectra with anal- 
ogous results for a sample of mock Lyman a forests. 

We used simulations run with t he parallel hydrody- 
namical (Tree SPH) code GADGET-2 JSpringel et al.l l200ll: 
ISpringell [20051. The simulations were performed with pe- 
riodic boundary conditions with an equal number of dark 
matter and gas particles and used the conservative 'entropy- 
formu lation' of SPH proposed by ISpringel fc Hernauisu 
(2002). Radiative cooling and heating processes were fol- 
lowed for a primordial mix of hydrogen and helium. We as- 
sumed a mean UV background produ c ed by quasars and 
galaxies as given by lHaardt fc Madaul Jl996T) with helium 
heating rates multiplied by a factor 3.3 in order to fit ob- 
servational constraints on the tem perature evolution of the 
IGM. More details can be found in lViel et ail j2004alL 

The cosmological model corresponds to a 'fiducial' 
ACDM Universe with parameters Qom = 0.26, SIoa = 
0.74, f2 g b = 0.0463 and Hp = 72 km s _1 Mpc -1 (the B2 
scries of lViel et aLll2004af) . We have used 2 x 400 3 dark mat- 
ter and gas particles in a 120 /t _1 comoving Mpc box. The 
gravitational softening was set to 5 hT 1 kpc in comoving 
units for all particles. We note that the parameters cho- 
sen here, including the thermal history of the IGM, are in 
perfect agreement with observational constraints including 
recent results on the CMB and other results obtained by 
the Lyman a forest community (e .g. ISpergel et alj 120031 : 
IViel et alj|2004atlSeliak et alJl2005l) . 

The z = 1.8 output of the simulated box was pierced to 
obtain 50 triplets of lines of sight carefully reproducing the 
observed triplet mutual separations. The same was done for 
50 sestets of lines of sight reproducing the observed sestet 
and 50 pairs of lines of sight at the same angular separation 
as Pair U. Fifty different realisations of Pair A spectra and 
50 of Pair Q were obtained from the output box at redshift 
2 = 2.4. 

The velocity extent of the spectra is of 10882.3 km s _1 
for the sestet, the triplet and Pair U, and of 11683.9 km s _1 
for Pair A and Pair Q. The pixel velocity size is fixed for each 
spectrum and slightly increases with the redshift of the box, 
being dv ~ 2.66, 2.85 km s _1 for z = 1.8, 2.4 respectively. 
These values are comparable with the average pixel size of 
the observed spectra, du bs — 2.64 km s _1 . 

The mean H I optical depth in the simulation box is de- 
termined by the UV background at the considered redshift. 
The corresponding observed quantity is the so called effec- 
tive optical depth, r c g, which is obtained from the average 
flux measured in QSO absorption spectra: / = exp(— r c g). 
The T e ff of the simulated spectra can be rescaled to a given 
value simply by multiplyi ng the optical depth of every pixel 
by a constant factor (e.g. lTheuns et alJll998bl:lBolton et al. 
2005). So for the spectra reproducing our sample of QSOs 
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Table 3. Meaning of the symbols used in the following figures 





LP sample pair sample 


simulations 


Auto correl. 


cross solid triangle 


empty triangle 


Cross correl. 


solid square 


empty square 


Cross corr. coeff. 


solid dot 


empty dot 



we have adopted the average values measured in the ob- 
served spectra: T e s = 0.10 for the triplet and Pair U (av- 
erage Lyman a redshift ZLya — 1.83), 0.12 for the sestet 
(z L ya ^ 1.83), and 0.21 for Pair A and Pair Q (z Lya ~ 2.34). 
The chosen values are within 1 per cent of the original values 
of the simulated spectra. In the computation of the corre- 
lation functions, we have used the average fluxes obtained 
as the mean of the flux values along the line of sight as for 
the observed spectra. We checked that the difference in the 
resulting correlation functions taking the true values for the 
average flux is less than 5 per cent. 

Finally, we added to the simulated spectra a Gaussian 
noise in order to reproduce the observed average S/N (per 
pixel): S/N — 5 for the triplet, the sestet and Pair A; S/N — 
9 for Pair U and S/N = 15 for Pair Q. 



3 CORRELATION PROPERTIES OF THE 
IGM: GENERAL PROCEDURE 

On the basis of the interpretation of the Lyman a forest as 
due to a continuous density field with a univocal correspon- 
dence between density and transmitted flux, we computed 
the correlation properties of the transmitted flux in QSO 
lines of sight and regarded them as indicators of the corre- 
lation properties of matter in the IGM. 

We selected in each normalised spectrum the region be- 
tween the Lyman /3 emission (or the shortest observed wave- 
length, when the Lyman /3 was not included in the spectrum) 
and 5000 km s" 1 from the Lyman a emission (to avoid prox- 
imity effect due to the QSO). Absorption lines due to ions 
of elements heavier than hydrogen 'contaminate' the Lyman 
a forest and c an give spurious contribution to the cluster- 
ing signal fsee lKim et al.ll2004l for a discussion in the case 
of single lines of sight). We flagged and removed the spec- 
tral regions where metal lines and Lyman a absorptions of 
damped and sub-damped systems occurred inside the Ly- 
man a forest. 

Given the normalised transmitted flux, /, as a func- 
tion of the velocity v» along the line of sight and the an- 
gular position 9 on the sky, we define 5f = (/ — /), where 
the average flux, /, is computed for every spectrum as the 
mean of the transmitted flux over all the considered pix- 
els in that spectrum. We neglected the redshift evolution of 
the average transmitted flux in the Lyman a forest of the 
individual spectra, which translates into the re dshift evo- 
lution of the mean Hi opacity of the Un iverse jKim et alJ 
120021: ISchave et alJl20oA IViel et al.ll2004aTl . because we ver- 
ified that its effect on the correlation function is negligible. 

Table [3] provides details about the symbol convention 
used to plot our results. 




100 1000 
Av. (km s" 1 ) 



Figure 3. Comparison of the auto correlation function for our 
sample of observed QSO spectra (solid triangles), for the LP sam- 
ple of observed QSO spectra (crosses, shifted in velocity for clar- 
ity) and for the simulated sample of spectra (empty triangles) as 
a function of velocity separation along the line of sight. Error bars 
on are highly correlated and they are plotted only to give an 
idea of their extent. 

4 THE FLUX AUTO CORRELATION 
FUNCTION 

The unnormalised auto correlation function of the flux along 
the line of sight is defined as: 

^(Av l{ ) = (5 f (v h e)5 f (v l} + Av n ,e)), (3) 

following previous st u dies on the same subject (e.g. 
[ McDonald et alJ l2000t iRollinde et"aH 120031: iBecker et all 
I2004I) . The auto correlation function for our sample of QSO 
spectra was obtained by averaging over all the pixels of all 
the QSOs. The results were binned in 50 km s _1 velocity 
bins. The 1 a error on this measure is estimated extracting 
50 samples of 15 objects from our sample with a bootstrap 
method, computing the auto correlation function for each 
sample and determining the standard deviation of the dis- 
tribution. 

The obtained correlation function is in v ery good agree- 
ment with the analogous determination by IRollinde et alJ 
(2003) and in qualitative agreement with the Lyman a 
lines two point correlation function at similar redshifts 
<Kim et al.ll200 ll) . 

The same procedure was adopted to compute £jj for the 
QSO spectra of the control sample taken from the UVES LP. 
These spectra have a much larger S/N (~ 50 in the Lyman 
a forest region) than our spectra, they have been selected 
to be free from damped Lyman a absorptions and they do 
not belong to known close QSO groups. 

The auto correlation function for the simulated spec- 
tra was computed as the arithmetic mean of the correlation 
functions obtained for 50 realisations of the observed sample 
and the error is the corresponding standard deviation. It is 
important to recognise that the computed error bars both 
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for the observed and simulated £^ are strongly correlated. 
This is due to the fact that every pixel contributes to the 
correlation function in several velocity bins. 

In Fig. [3]we show the comparison of the three auto cor- 
relation functions in velocity space, error bars are included 
just to give a qualitative idea of their size. 

The agreement at the 3 a level is good. However, we 
would like to point out two kinds of discrepancies at the 
1 a level which, in our opinion, are the signatures of two 
interesting effects. 

1. Both observed correlation functions have an amplitude 
in the first bin (centred at 50 km s _1 ) which is larger than 
the simulated one. We ascribe this to a scantiness of strong 
lines (log ./V(HI) > 15.5) in the simulated spectra (see e.g. 
iTheuns et aljEoOSi ) due to the limited extent of simulations 
in comoving space, which implies that large over densities 
(giving rise to strong Lyman a lines) are not appropriately 
probed. We found support to this hypothesis from the anal- 
ysis in Sections 4.1 and 6.1.1. 

2. The correlation function computed with our sample of 
QSOs tends to show a stronger clustering than the simu- 
lated one up to large separations, while the LP correlation 
function is compatible with the mock correlation function 
already at the second bin. The main difference between the 
two samples is the fact that our QSOs were selected to be- 
long to close pairs and associations, thus they could trace 
over-clustered regions. Although we excluded the spectral 
portion within 5000 km s" 1 of the QSO emission, an en- 
hanced clustering signal could still be present. We discuss 
this hypothesis in Section 4.2. 

We notice that cosmic variance could also be at the ori- 
gin of these small differences. Indeed auto correlation func- 
tions of single observed QSO spectra show large variations 
among them and we are working with relatively small sam- 
ples. On the other hand, all the simulated spectra have been 
drawn from the same box of 120 h -1 comoving Mpc, which 
although notably large for a si mulation, possibly unde resti- 
mates the cosmic variance (see IMcDonald et al.ll2000t for 
possible ways of correcting the simulated correlation func- 
tion for errors induced by the finite box size). 

4.1 Contribution of the stronger absorptions 

Lyman a lines of increasing Hi column density correspond 
to structures of increasing overdensity and smaller charac- 
teristic size. Going from log AT (Hi) ~ 12 to logiV(Hl) > 15 
we probe from under-dense, modestly over-dense regions, to 
highly over-dense filaments and their intersections. We ex- 
pect that the major contributors to the clustering signal are 
the stronger absorption lines, which trace the higher den- 
sity peaks. The aim of this section is twofold. On the one 
hand, by excluding the stronger lines we want to verify the 
clustering properties of the gas in the 'true' IGM, far from 
the densest regions. On the other hand, we want to test 
the relative behaviour of the observed and simulated auto 
correlation functions without those lines. 

Note that a careful study of the impact of strong ab- 
so rption systems on t he fl ux power spectrum was performed 
bv lViel et all (l2004bl) and IMcDonald et al.1 J2005I) . 

To select the strong absorption lines without carrying 
on the time consuming Voigt fitting procedure of the ob- 
served and simulated spectra, we can relate the transmit- 




Av, (km s _1 ) 

Figure 4. Same as Fig.[!|]but selecting only the pixels with flux 
value / > 0.1 both in the observed and simulated spectra. This is 
equivalent to exclude the lines with column density log N(H i) > 
13.96. 



ted flux values to a corresponding column density assum- 
ing a typical Doppler parameter. We eliminated the pix- 
els with / < /thres = 0.1 or thi > 2.3. This approxi- 
mately corresponds to excluding the absorption lines with 
logiV(Hl) > 13.96, if we adopt the formula for the optical 
depth at the centre of the line: 



T0,Ly-a 



7.58295 x 10 



-is N(c 



b (km s 



(4) 



and a Doppler parameter b — 30 km s 1 chara cteristic of 
the L yman a forest absorptions at redshift z ~ 2 ijKim et alJ 
1200 ID . 

This technique is efficient for data with S/N > 
l//thrcs = 10, so we applied the selection only to the LP 
QSO sample and to the sample of simulated spectra for 
which, in this case, we adopted a Gaussian noise giving 
S/N = 50. 

Fig. 2] shows the resulting correlation functions. The 
amplitude of the observed £jf has reduced to about one third 
of the initial value but still shows a significant clustering 
signal. This is an indication that matter is still clustered 
also far from the most over-dense regions. The threshold on 
thi, transformed with eq. (2), corresponds to p/p < 6.5. 

What is worth noting is that the amplitude of the sim- 
ulated correlation function has decreased of about one half 
and is now fully compatible with the observed correlation 
function. This suggests that strong lines are deficient in the 
simulated spectra so that their contribution to the correla- 
tion function is smaller than for the observed spectra. When 
they are excluded from the computation the two correlation 
functions become consistent. 
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4.2 Selection effects in the observed spectra 

There is the possibility that our observed sample is affected 
by a selection bias since the QSOs were required to be at 
close angular separations and, in most cases, they have sim- 
ilar emission redshifts. As a consequence, the observed lines 
of sight could be biased toward more clustered regions than 
average. 

We tried to estimate this effect on the auto correla- 
tion function by excluding larger and larger portions of 
the observed spectra from the QSO emissions (up to 20000 
km s _1 ). No significant decrease of the amplitude of the 
correlation function has been detected. 

This test could prove inadequate as far as the ses- 
tet is concerned since the redshifts of some member QSOs 
are separated by up to ~ 50000 km s" 1 (see Fig. 0. 
As a consequence, our measure c ould be affected by 
the transverse proximity effect (e.g. Ijakobsen et alJ |2003t 
IWorseck fc Wisotzkill2006h or by possible clustering associ- 
ated with the presence of a foreground QSO. Indeed, in the 
case of Pair Q, we found a damped Lyman a system in the 
higher redshift QSO spectru m at the redshift of the other 
QSO JP'Odorico et al.ll2002l) . From a preliminary analysis 
of the sestet, strong Lyman a absorptions (log iV(Hl) ~ 15) 
are observed in the spectra of S3, S5 and S6 at ~ 600 — 750 
km s _1 from the emission redshift of the lower redshift QSO, 
SI, and a strong Lyman a with associated Civ absorption 
is present in the spectrum of S2, the furthest away QSO, 
at ~ 1000 km s" 1 from the redshift of SI. Considering the 
other foreground object, S4, Lyman a absorptions within 
200 km s _1 of its redshift are observed along the line of 
sight to S2 (log AT (Hi) ~ 14) and S3 (log AT (Hi) ~ 15) with- 
out associated metals detected. As already mentioned in the 
Introduction, we will carry out a detailed study of the co- 
incidences among multiple lines of sight in a future paper, 
with data of improved S/N. 

In the following, we will use the auto correlation func- 
tion determined from the LP sample to avoid systematic 
uncertainties introduced by the above mentioned possible 
selection effect. 



5 CLUSTERING OF THE FLUX ACROSS THE 
LINE OF SIGHT 

In this section we exploit the potentialities of our sample of 
QSO pairs by determining the clustering properties of the 
IGM across the lines of sight. The great advantage with re- 
spect to the correlation function along the line of sight, in 
particular for a sample like ours showing a large variety of 
pair separations, is that we have the guarantee of sampling 
true spatial separations between the pixels, the effect of pe- 
culiar velocities being negligible or absent. 

5.1 The cross correlation function in redshift 
space 

As a first approach, we computed the cross correlation func- 
tion extending in a natural way the procedure adopted for 
the auto correlation function. 

Every pixel along the line of sight is considered as an 
element of the density field at the QSO angular position in 



the sky and at a distance from the observer (comoving along 
the line of sight) corresponding to the redshift of the pixel: 

where E(z) is defined as: 



E{z) = v /n 0m (l + z) 3 +OoA. (6) 

In eq. (5), there is the implicit hypothesis that peculiar 
velocities contribute negligibly to the measured redshift in 
the Lyman a forest. This state ment is supported b y recent 
measurements using QSO pairs llRauch et al.ll2005T) . and we 
will see a posteriori that the comparison of the auto and 
cross correlation functions makes this procedure allowable. 

The cross correlation function of the transmitted flux 
between two lines of sight at angular separation A9 is defined 
as: 

££(Ar) = (8 f {0,r hl )8f(B + A9,r h2 )), (7) 

where, Ar — \fiy?, 1 + r| 2 — % r \\,i r \\,2 cos A8) is the spatial 
separation between pixel 1 at 7*11 x along one line of sight and 
pixel 2 at rii 2 along the paired line of sight. 

The error on the correlation function cannot be esti- 
mated with a bootstrap technique, as done before, due to 
the limited number of pairs contributing to each separation 
(in particular at the smaller separations). In order to give 
an evaluation of the significance of the observed signal, the 
same cross correlation function was computed replacing one 
QSO in each pair with a control QSO, then repeating the op- 
eration replacing the other QSO of every pair with another 
control QSO. The control sample is formed by the LP QSO 
spectra. We chose LP QSOs at redshifts close to those of the 
original objects and shifted the spectra in order to match it 
exactly. We derived an indicative error as the r.m.s. devia- 
tion from zero of the two control functions in two regions: 
Ar < 5 h~ x Mpc, for which the uncertainty is larger, and 
Ar > 5 h' 1 Mpc. 

We computed the cross correlation function also for the 
sample of mock spectra. The simulated spectra are char- 
acterised by the redshift of the output box and a velocity 
extent (see Section 2.2). In order to assign a redshift value to 
every pixel, we gave the central pixel of every spectrum the 
redshift of the corresponding output box, then we numbered 
the pixels one by one transforming the velocity size of the 
pixel into a redshift size. Once the redshifts were determined 
pixel by pixel, we followed the same procedure adopted for 
the observed spectra for the 50 simulated samples and com- 
puted the average cross correlation function and its 1 a stan- 
dard deviation. 

The result of our computation is shown in Fig. com- 
pared with the observed cross correlation function of the 
pairs and of the control sample. There is a very good agree- 
ment between the two functions. Only in the first bin there 
is an indication of a lower value for the simulated cross cor- 
relation with respect to the observed one. This is consistent 
with the result obtained for the auto correlation function if 
we take into account the fact that the first bin of corre- 
sponds to the second bin of fjf . 
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Figure 5. Cross correlation function of the transmitted flux for 
our sample of pairs and groups of QSOs (solid squares) com- 
pared with the analogous correlation function computed with the 
sample of mock spectra (empty squares). The error bars on the 
simulated data are 1 a standard deviations for 50 simulated sam- 
ples while the error bars on the observed data have been obtained 
from the r.m.s. deviation of the two control samples plotted as 
dashed lines (see text). Simulated data arc slightly shifted in Ar 
for clarity. 

5.2 The cross correlation coefficient 

A measure of the transverse clustering properties of the 1GM 
which is less affected by peculiar velocities is the flux cross 
correlation coefficient, 

x{(Ae) = (5f(6,v }l )8 f (9 + A9,v ll )), (8) 

where, every pixel along one line of sight is correlated with 
the one face-to-face in redshift space along the paired line 
of sight and the result is averaged over all the pixels in the 
common redshift interval. 

Every pair of QSOs at angular separation A9 gives one 
value of Xx (A9), and a sample with several pairs at different 
separations, as is our sample, gives an estimate of the corre- 
lation function. At a given redshift, the angular separation 
A6 corresponds to a velocity separation Av± = cF(z)A9, 
where c denotes the speed of light, and F(z) is a dimension- 
less function of redshift that includes all the dependence on 
the global cosmological metric, tn the cosmological model 
that we have adopted: 

E(z) r[dz/E(z)} 
F(z) = '\ (9) 

and E(z) is defined in eq. (6). 

We computed the range of velocity separations covered 
by each of our pair of spectra then we grouped the pairs in 
velocity bins of variable width and computed the average 
cross correlation coefficient for every group. Given the small 
number of pairs in every group (a maximum of 3 QSO pairs) 
the uncertainties reported on these determinations are com- 
puted by applying a simple error propagation to eq. (8), 



Figure 6. Comparison of the cross correlation coefficients for 
our sample of observed spectra (solid dots) and of simulated ones 
(empty dots) as a function of the velocity separation, Av± = 
cF(z)A9, corresponding to the angular separation, A9, of the 
QSO pairs. The angular separation computed at z = 1.8 is re- 
ported in the top axes. Observed values are shifted by 10 km s -1 
for clarity. The lower observed value at Av± ~ 500 km s~ 1 was 
obtained excluding pair S1S3 and removing a strong coincident 
Lyman a line in the T1T3 pair (see text). Error bars on the ob- 
served values along the x axis represent the velocity range covered 
by the considered pairs. Errors on the \x represent the uncer- 
tainty on that measure for the observed values and the cosmic 
variance for the simulated values (see text). 

thus they account only for the pixel statistics and the noise 
associated with the transmitted flux but they are not repre- 
sentative of the true error due to the cosmic variance. 

In the case of the simulated spectra, we had 50 realisa- 
tions of each of our QSO pairs so we could obtain in every 
velocity interval defined for the observed pairs an average 
cross correlation coefficient with its error, that in this case 
is the standard deviation of the distribution of values. 

Results are shown in Fig.|S] The two points at the small- 
est separations are given by Pair U + S5S6 and Pair A, 
respectively. The two samples of values are in very good 
agreement except for the bin at Av± ~ 500 km s _1 (~ 5 ar- 
cmin at z ~ 1.8) for which the observed value is ~ 7 a above 
the simulated value. The observed Xx (AO) in this velocity 
bin is obtained by averaging two pairs: S1S3 and T1T3 for 
which Xx - 0.040 ± 0.003 and 0.0177 ± 0.0015, respectively. 

The two pairs of spectra were searched for peculiar fea- 
tures that could boost the signal by selecting those pixels for 
which 5f(9,v\\)5f(6 + A9,v\[) > 0.4. fn Fig. Qwe show the 
results of this research. In the case of SI (up left plot) about 
44 per cent of selected pixels have values / > 2 or / < 0, 
that is they are strongly affected by noise (in the case of S3 
only ~ 14 per cent of pixels are in those range of values). As 
a consequence, the large value of the cross correlation coeffi- 
cient for pair S1S3 could be due to the low S/N in the spec- 
trum of SI and cannot be determined reliably. On the other 
hand, in the case of pair T1T3, Fig. □ (centre plots) shows 
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Figure 7. Selection of the pixels giving values Sf(9, v^)Sf(0 + 
A0,U||) > 0.4 for the pairs S1S3 (left), T1T3 (centre) and 
T2T3 (right). The two upper plots show the values of the flux 
for the selected pixels in the indicated QSO while the lower plot 
report the value of the above product. 




Ar (comoving Mpc) 

Figure 8. Cross correlation function of the transmitted flux, , 
for our sample of pairs and groups of QSOs (squares, slightly 
shifted in Ar) compared with the observed cross correlation coef- 
ficients, Xx (dots), as a function of spatial separation (see text). 
Error bars both on an d nave been determined from simu- 
lations. 



a feature giving a strong signal at A ~ 3300 A. We identi- 
fied in the spectra of Tl and T3 two coincident strong H I 
Lyman a absorptions (logiV(Hl) ~ 15) at this wavelength, 
the former with a clear associated C IV doublet while for the 
second there is the possible detection of the C ivA 1548 line 
only. No corresponding H I Lyman a line is observed at this 
redshift along the third line of sight of the triplet, T2. If 
we mask the redshift interval covered by the coincident lines 
and recompute x{ for T1T3 we get: 0.0075 ± 0.0015, which 
is in better agreement with the simulated value as shown in 
Fig. 

In the right hand plots of Fig. [7| we show the case of 
pair T2T3 (at 8.9' angular separation) where a coincident 
absorption system is present with characteristics very simi- 
lar to the one in pair T1T3, and the cross correlation coeffi- 
cient has a value consistent with zero. This suggests that the 
presence of a coincident absorption system could be a neces- 
sary but not sufficient condition to explain a large value of 
Xx • It is interesting to note that also the recently published 
two point correlat ion function of the C IV a bsorptions in the 
LP QSO spectra dScannapieco et ai1l2006t) shows an excess 
of clustering signal at ~ 500 km s _1 whose origin has not 
been satisfactorily explained. 

New observations are going to be carried out with 
UVES to increase the S/N in the two QSO pair spectra 
S1S3 and T1T3, in order to confirm the large value of the 
cross correlation coefficient and possibly detect other coinci- 
dent C IV systems. On the other hand, we will try and collect 
other QSO pairs at the same velocity separation to increase 
the statistics and we will look into physical mechanisms that 
could explain our results. 

In Fig. El we compare the cross correlation coefficients 
with the cross correlation function. The angular separation 



A6 between two QSO lines of sight was transformed into a 
comoving spatial separation, Ar, with the formula: 

cA9 dz' , , 

where E(z) was defined in eq. (6). Considering the error 
bars computed from the simulated pairs, there is a good 
agreement within 1 a. The large variations from one data 
point to the other in the cross correlation coefficients should 
be due mainly to the small number of QSO pairs (between 
1 and 3) contributing to each point. On the other hand, the 
smoothness of f£ is artificially increased by the fact that the 
values in the different bins are not independent. 



6 DISCUSSION 

From the cosmological point of view, one of the most in- 
teresting and challenging applications of the kind of cal- 
culations carried out in this paper is the constraint on 
the geometry of the Universe (in particular, the esti- 
mate of 57qa h~ 2 ) by linking angular separations and red- 
shift differences in the hypothe sis that the observed cor- 
relation properties are isotropic (|Alrock_fc^aczynski||l97^; 
iMcDonald fc Miralda-Escudelll999l iHui. Stebbins fc Burlesl 



11999 

A large number of QSO pairs at different angular sepa- 
rations ( see lMcDonaidi r2003 ) is needed to obtain a measure 
of the cosmological parameters as accurate as those recently 
produced by the WMAP team dSpergel et alJl2003ft . This is 
why we do not attempt to derive a measure of fioA with the 
present sample. 

Nevertheless, a comparison of the observed auto and 
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Table 4. Effective optical depth adopted to build the three sam- 
ples of mock spectra used to compute the correlation functions in 
Figs [Till and HT1 



Ar (comoving Mpc) 

Figure 9. Comparison of the cross correlation function for our 
sample of QSO pairs (squares) with the auto correlation func- 
tion computed for the LP QSO sample (crosses) as a function of 
comoving spatial separation across and along the line of sight, 
respectively. The cross correlation function is slightly shifted in 
Ar for clarity 



cross correlation functions is interesting to qualitatively 
evaluate the distortion due to peculiar velocities and the 
correctness of the adopted cosmological parameters. The 
two functions are shown in Fig. [5] Note that, due to the 
minimum angular separation between our QSO pairs, there 
is no transverse clustering signal corresponding to the first 
bin of the auto correlation (centred at Av = 50 km s _1 or 
Ar ~ 740 kpc). The correspondence between £^ and ^ is 
very good suggesting that the adopted cosmological param- 
eters are reasonable and that peculiar velocities do not play 
a major role in the IGM gas on velocity scales > 100 km s _1 . 

In the next section, we survey the behaviour of the sim- 
ulated correlation functions when varying three parameters 
characterising the physical state of the gas. 



6.1 Exploring the parameter space of the 
simulations 

6.1.1 Effective optical depth 

As already explained in Section 2.2, the effective optical 
depth of a simulated spectrum can be rescaled to a new 
value simply by increasing or decreasing the optical depth 
of every pixel in the spectrum by a constant factor, which is 
equivalent to increase or decrease the column density of the 

absorption lines in the spectrum. 

Considering previous obse rvations (|Vie l et a I] l2004al) 
and the relations computed bv lSchave et alj " (l200^) for the 
redshift evolution of the effective optical depth, we derived 
3 a upper and lower limits for r e ff at the redshifts of the sim- 
ulated spectra, which are reported in Table[I] Then, we built 
two new samples of mock spectra starting from the original 
spectra produced from the simulation, one with the mini- 
mum and one with the maximum values of r e ff and recom- 
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Figure 10. Auto correlation function for the LP sample (crosses) 
related to the average auto correlation functions of three sam- 
ples of simulated spectra with different T e ff . The long dashed line 
traces the correlation function for the sample with the minimum 
optical depths, the dashed line is for the sample with the ob- 
served optical depths and the solid line is for the sample with the 
maximum optical depths (see Table HI . 



puted the auto and cross correlation functions. Each sample 
is formed by 50 realisations of the observed sample of QSO 
pairs and we determined the average value for the correla- 
tion function. 

Our results are shown in Figs llOl and llll At the small- 
est scale, probed only by the auto correlation function, the 
agreement between observed and simulated clustering signal 
improves adopting the larger effective optical depths. How- 
ever, at the Jeans scale, corresponding to the first bin of f£ 
and to the bin centred at 100 km s _1 of the larger t c s 
possibly overestimates the clustering amplitude. 

Since the effective optical depth measured in the orig- 
inal simulation is in very good agreement with the one ob- 
servationally determined from QSO Lyman a forests (see 
Secion. 2.2), we infer that the discrepancies between simu- 
lations and observations are not due to a wrong estimate of 
Tea but, as already mentioned in Section 4, to the deficiency 
of strong lines in the simulations. In the spectra with the 
larger r c ft the column density of all the lines is increased 
and the average flux is decreased, causing a general increase 
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Figure 11. Cross correlation function for our sample of QSOs 
(solid squares) compared with the average cross correlation func- 
tions of three samples of simulated spectra with different r e ff. 
The long dashed line traces the correlation function for the sam- 
ple with the minimum optical depths, the dashed line is for the 
sample with the observed optical depths and the solid line is for 
the sample with the maximum optical depths (see Table HI 

of the amplitude of the correlation function, partially com- 
pensating the small number of strong lines. 

In order to verify our statement, a quantitative analysis 
of the difference in the Lyman a column density distribution 
function (the number of absorption lines per unit column 
density per unit redshift) of obs erved and simulated spectra 
(as done bv lTheuns et alJl2002l for small simulated boxes) 
is necessary. Furthermore, more pairs at small separations 
(< 1') would allow to decrease the error bars and test if the 
behaviour of clustering is the same as along the line of sight. 

6.1.2 Temperature and 7 

The influence on the flux cross correlation function of the 
temperature and of the power- law index of the temperature- 
density relation adopted for the gas (see eq. (1)) are ex- 
plored. The aim of this analysis is to point out the modifi- 
cations of the simulated correlation function when the pa- 
rameters are varied, and not to derive precise constraints on 
the same parameters from the comparison with observations. 
With a running time of 2-3 weeks for simulation, it is impos- 
sible to run an extensive parameter study for different ther- 
mal histories, thereby we decided to rescale a posteriori the 
temperature-density relation at a given output by assuming 
that all the gas particles obey this relation. We then recom- 
puted the ionization fractions for each gas particle with the 
proper UV background. This approximation does not ac- 
count for the effects that the corresponding change in the 
gas pressure would have on the gas distribution. However 
we explicitly checked that the differences in the correlation 
function for the fiducial run and for a rescaled temperature 
relation with the same T — 7 are of the order of f per cent 



Figure 12. Cross correlation function for our sample of pairs 
(solid squares) compared with the average cross correlation func- 
tions of samples of simulated spectra with different temperature 
and exponent 7 of the temperature-density relation (eq. (1)). 
Short dashed lines are for simulations in which 7 = 1.6 and tem- 
perature is T = 8000, 22000 and 50000 K going from the upper 
to the lower curve. The long dashed lines represents simulations 
with 7 = 1.1 while the temperatures are the same as before but 
they increase from the lower to the upper curve 



(see IViel et aljl2004al . for the effect of the rescaling on the 
flux power spectrum). 

We built six sets of simulated spectra each one formed 
by 50 realisations of the observed sample of pairs. Every set 
is characterised by a value of 7 (1.1 or 1.6) and a value of the 
temperature (8000, 22000, and 50000 K), while the effective 
optical depth is the one used for the main sample (refer to 
Section 2.2). 

In Fig. 1121 we plot together the observed f£ with the av- 
erage correlation functions computed from the six groups of 
simulated spectra. At the smallest separation probed by the 
observations, the predictions of the different models are well 
separated and could be tested with a larger sample of close 
QSO pairs. The models with 7 = 1.1 give a larger cluster- 
ing signal than those with 7 = 1.6 for all temperatures and 
are in better agreement with the data at the small scales. 
An improved concordance between observations and simu- 
lations for lower values of 7 was observed also in the case 
of the flux probability distribution function estimated from 
the Lyman a forest (James Bolton, private communication). 

The exponent of the temperature-densit y relation is de- 
termin ed by the ionization history of the gas. lHui fc Gnedlnl 
il997fl demonstrated that instant reionizations of H I occur- 
ring at redshifts decreasing from z — 10 to 5 would imply 
shallower and shallower 7 at redshifts z ~ 3. However, also 
for very late reioniz ations 7 should have increased a bove val- 
ues ~ 1.4 by 2 ~ 2. iRicotti. Gnedin fc Shuil fcOOd) studied 
the evolution with redshift of the equation of state of the gas 
and showed that a reionization of He 11 occurring at z ~ 3 
would cause a sudden increase in the gas temperature and 
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a corresponding decrease in the value of 7 (their fig. 13). 
Observationally, this has been studied mainly by using the 
distribution of Doppler para meters and column densities of 
the L yman a forest lines fe.g. lSchave et alj |2000; Kim e^y 
| 2002j). but also wit h metal absorptions fe.g Tsongailall^9^ 
IVladilo et alJ l2003). Our result supports previous observa- 
tional evidences of a second reheating of the Universe hap- 
pening at z ~ 3, most likely due to the reionization of Hell 
and it requires confirmation by a larger sample of observed 
spectra and a set of more refined mock spectra. 



7 CONCLUSIONS 

In this paper, we exploited the capabilities of the largest 
sample of high resolution UVES spectra of QSO pairs to 
study the 3-dimensional distribution properties of baryonic 
matter in the IGM as traced by the transmitted flux in the 
QSO Hi Lyman a forests. Our sample is formed by 21 QSO 
pairs evenly distributed between angular separations of ~ 1 
and 14 arcmin, with Lyman a forests at a median redshift 
z ~ 1.8. We selected also 8 UVES QSO spectra from the 
ESQ Large Program 'The Cosmic Evolution of the IGM' 
toergeron et alJ|2004h to compute the correlation function 
along the line of sight and to be used as a control sam- 
ple for the cross correlation function (see Section 5.1). We 
compared the observed sample with a set of mock spectra 
drawn from a cosmological hydro-simulation run in a box 
of 120 h" 1 comoving Mpc, adopting the cosmological pa- 
rameters of the concordance model. The simulated sample 
reproduces 50 different realisations of the observed sample. 
In the following, we resume our main results: 

1. The clustering properties of matter in the IGM are the 
same in the direction parallel and transverse to the line of 
sight when using the parameters of the concordance cosmol- 
ogy to transform the angular distance into velocity separa- 
tion. As an implication, peculiar velocities in the absorbing 
gas are likely smaller than ~ 100 km s _1 . 

2. Matter in the IGM is clustered on scales smaller than 
~ 200 km s _1 or about 3 h~ J comoving Mpc. We verified 
that the clustering signal is significant also for the slightly 
over-dense gas (thi > 2.3 or p/p < 6.5) although on smaller 
scales. 

3. The simulated correlation functions are consistent with 
the observed analogous quantities at the 3 a level, although 
they systematically predict lower clustering at the smaller 
scales. The agreement becomes better when only the lower 
density regions are selected for the computation or when the 
effective optical depth of the simulated spectra is fixed to a 
larger value (marginally consistent with previous extensive 
observational results on the redshift evolution of the effective 
optical depth) . These are indications of a deficiency of strong 
absorption lines in the simulated spectra that needs further 
investigation. 

4. We observed an improved consistency between observa- 
tions and simulations also when a lower 7 is adopted in the 
equation of state of the gas, T = To (p/p) 7_1 (7 = 1-1 in- 
stead of the standard 1.6). This result hints to a late Hell 
reionization epoch whose effects on the IGM could still be 
measured at the redshifts investigated here. 

5. We measured an enhanced clustering signal for the cross 
correlation coefficient at a transverse velocity separation 



A?;x ~ 500 km s _1 . This velocity distance matches the line 
split of the Civ doublet. We propose to gather more obser- 
vational and theoretical material to shed light on this result. 

From the analysis carried out in the present study, we 
evidenced the need of increasing our sample of observed QSO 
pairs, in particular at the small angular separations (~ 1 ar- 
cmin), but also at Av± ~ 500 km s _1 in order to confirm the 
overdensity of coincident Lyman a lines at this transverse 
distance. On the other hand, we realised the necessity of a 
dedicated effort on the comparison between observed and 
simulated QSO spectra, in particular as far as the statistics 
of absorption lines is concerned (number density, Doppler 
parameter and column density distribution, etc.). 
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